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ABSTRACT 

We define the notion of "diffusion algebras". They are quadratic Poincare- 
Birkhoff-Witt (PBW) algebras which are useful in order to find exact expressions 
for the probability distributions of stationary states appearing in one-dimensional 
stochastic processes with exclusion. One considers processes in which one has N 
species, the number of particles of each species being conserved. All diffusion al- 
gebras are obtained. The known examples already used in applications are special 
cases in our classification. To help the reader interested in physical problems, the 
cases N = 3 and 4 are listed separately. 
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1 Introduction. 



One-dimensional stochastic processes with random-sequential updating have sta- 
tionary states (far away form equilibrium) with very interesting physical properties. 
One observes phase transitions about which still little is known. As opposed to 
equilibrium states where the phase transitions are essentially a bulk phenomenon, 
in stationary states, the boundary conditions play an essential role [jl| . One observes 
bulk induced phase transitions if one looks to problems on a ring [^, boundary in- 
duced phase transitions in the case of open systems , or a combination of both 
[^]. When dealing with phase transitions, one is interested in exact expressions for 
the relevant physical quantities like the current densities (which are zero for equi- 
librium problems) and various correlation functions. A major step in this direction 
was achieved when it was understood that in certain cases one can use the so-called 
"matrix product states" approach [§, |l|]. In our opinion, from a mathematical point 
of view, this approach was loosely defined (for an illustration see Appendix A). The 
aim of this paper is to improve upon this situation and, as a bonus, to give a large 
number of processes where one can find " matrix product states" . We define a class 
of quadratic algebras to which we have coined the name "diffusion algebras". We 
also show how to construct all of them. We believe that in this way one can bring 
some mathematical "beauty" in what was not, up to now, a systematic approach. 

We start with the physical problem in order to explain the motivation of our 
mathematical work. 

In the present paper we consider a restricted class of processes in which we 
take N species of particles with — 1 conservation laws in the bulk. We take a 
one-dimensional lattice with L sites and assume that one has only nearest-neighbor 
interaction with exclusion (there can be only one particle on a given site). In the 
time interval dt only the following processes are allowed in the bulk: 

a + p^p + a {a,p = 0,1,...,N -1) , (1.1) 

with the probability Qapdt {gap > 0). Eq.(|Ll|) is a symbolic equation used by 
physicists. Its meaning is that the particles a and (3 on successive sites exchange 
their places. Obviously the number of particles Ha of each species a are conserved 
(Z]«=o^ = L). One is interested in the probability distribution P(ai, ^2, • • • , «l), 
(ttfc = 0, 1, . . . , A^ — 1) for the stationary state. In order to obtain it, in the matrix 
product approach, one considers A^ matrices Va and A^ matrices acting in an 
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(1.2) 



[Da, Xp] 




In Eqs.( |1.2| ) {., .} represent anti-commutators. 

If one considers processes on a ring (periodic boundary condition), the un- 
normalized probability distribution has the following expression 0: 



Notice that the matrices Xa don't appear in Eq.( |1.3|) . The expression of the proba- 
bility distribution is special in at least two ways. As a consequence of the conserva- 
tion laws, Eq. ( |1.3|) connects only monomials with the same numbers no, ni, . . . , hn-i 
of generators Vo,Vi, . . . This implies that one can use different matrices for 

monomials for different values of the set hq, rii, . . . , n^v-i- This observation is impor- 
tant since for example, a given infinite-dimensional representation can have a finite 
trace for certain class of monomials but can diverge for another class of monomials. 
In order to use Eq.( |1.3| ), for the latter one can use a different representation which 
for example, is traceless for the first class of monomials but has a finite trace for the 
second class. This problem will be explained in detail in p. 

On the other hand if different representations can be used for all monomials 
one has the remarkable property that up to a factor (the expression ( |1.3| ) does not 
contain a normalization factor), the traces are independent on the representations 
one uses. For concrete calculations one takes therefore the representation with the 
smallest dimension. 

If one considers open systems, the bulk processes have to be completed by bound- 
ary processes (they break the conservation laws). On the left side of the chain (site 
1) and on the right side of the chain (site L) we assume that in the time interval dt, 
the particle a is replaced by the particle /5: 



P(ai, a2, . . . , ol) = Tr {Va^Va^ ■ ■ ■ V^^) 



(1.3) 



a 



/5, 



(1.4) 



with the probabilities 

L'idt , respectively RZdt . (1.5) 



The matrices appearing in ( |1.5| ) are intensity matrices for which the diagonal 
elements are given by the non-diagonal ones: 



= - ^ , = - ^ i?^ , (1.6) 



Py^a j3^a 
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The non-diagonal matrix elements are non-negative. 

For open systems, the un-normalized probability distribution functions are given 
by a matrix element in the auxiliary vector space: 

P(ai, a2, . . . , aL) = (0| V^J)^, . . .V^^ |0) , (1.7) 

where the bra (ket) state (0| (respectively |0)) are given by the following conditions: 

{Q\{LiVp + X^)=Q, {RiVp-X^)\Q) = Q, (1.8) 

The expressions ( |1.3| ) and ( |1.7| ) for the probability distributions are of httle use 
unless one finds matrices which satisfy the conditions ( |1.2|) and have a trace, or one 



finds matrices which satisfy ( [L.2D and have the property ( |1.8|) . In a very nice paper 
[IC] it was shown that for arbitrary bulk and boundary rates one can in principle 
construct infinite dimensional matrices which satisfy the conditions ( [L.2|) and (|1.8|) (if 
they have also finite traces is unclear). It is however practically impossible to obtain 
explicit expressions for these matrices. It is also not clear which supplementary 
relations besides ( p..2|) they satisfy. In other words, the algebraic structure behind 
the relations ( |1.2D is obscure. 

In the present paper we adopt a different approach. We start by searching for 
quadratic algebras which are of PBW type (this notion is explained in the beginning 
of Section 2) with a structure which is understood and look for those bulk rates for 
which a simplified version of the relations (|1.2| ) exist. Then we look to boundary 
conditions compatible with the simplified version of Eq.( p..8|) . We make the Ansatz: 



X„ = x„e, (1.9) 

where e acts clS cl unit element on 

eVa = V^e = Va , (1.10) 



and Xa are c-numbers. With this Ansatz, instead of Eq.( p..2| ) one obtains A^(iV — 1)/2 
relations for the matrices 

and instead of Eq. ( PT8| ) one obtains: 

{Q\{LiVp + x^e)=Q, {RiVp~x^e)\Q) = Q, (1.12) 
Taking into account the relations ( p..6| ), one has 

N-l 

5Z Xa = Q ■ 



0=0 
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The consequence of the Ansatz ( |1.9D will be that we will not be able to find 
expressions for the probability distributions for arbitrary bulk and boundary rates. 
Further limitations on the possible bulk rates will appear when we ask for the 
relations ( 1.11 ) to define quadratic algebras of PBW type. This implies that the 
ordered monomials 

Vl^^VlU-'-T^T-' (1-13) 

form a basis in the algebra. Algebras of this type will be called " diffusion algebras" . 

Once the algebras are known, we can ask if they are also useful for applications to 
stochastic processes. First one has to find for which algebras one can choose all the 



ga.13 non- negative. Next, since the relation ( |1 . 1 1| ) gives recurrence relations among 
traces of different monomials, one has to find in which cases one has not only the 
trivial solution for these recurrence relations (all the traces vanish). These cases can 
be used to study stochastic processes on a ring. Finally, one has to find for which 
boundary matrices one can find representations for which the conditions ( p..l2| ) are 
satisfied. This is not a trivial exercise. For the cases for which one finds solutions, 
one can get then the expression of the probability distribution for the open system. 

We would like to give a supplementary argument which motivated us to look for 
PBW algebras. Interesting physics appears when one has only infinite dimensional 
representations either for the ring problem or for the open system [||, 0, 0. The 
PBW algebras have at least one infinite-dimensional representation (the regular 
one), which doesn't mean that, in general, they don't have also finite-dimensional 
ones. 

In the present paper we concentrate on the diffusion algebras only. In a sequel 
[^] we will discuss in more detail the representation theory. We will also not look 
for physical applications except for stressing the cases when the rates Qafi can't be 
chosen non-negative. 

In Sec. 2 we will show that in order to have a diffusion algebra, the gap and Xa 
have to satisfy certain identities which are the equivalent of the Jacobi identities 
for the structure constants in the case of Lie algebras Q. It is remarkable that, as 
we are going to show in the next sections, we are able to find all the solutions of 
these identities. Inspecting the Eqs.( |l.llD one notices that if some of the 
not zero they can be rescaled in the definition of the generators T>a- The number 
of non-vanishing Xa will play an important role in the classification of diffusion 
algebras. 

In Sec. 3 we consider the case = 3 corresponding to the three species problem 



^Correctly speaking, we mean universal enveloping algebras of Lie algebras 
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(the case = 2 is well known |TT|). The case = 3 is not only interesting on its 
own but is relevant for understanding the case for arbitrary described in Sec. 4. 
We are going to recapture all the known examples p|, |l^, |l3l which are interesting 
for applications and get a few new algebras which are interesting on their own. 

We would like to mention that the case = 3 was approached previously from 
two other points of view. In Ref . , one has looked directly to the open system, and 
classified the type of boundary matrices appearing in Eq. (|1.12| ). This is possible since 
the boundary matrices are intensity matrices which, as explained in Ref. 0, have 
special properties. Next, one has looked for bulk rates (see Eq. (|l.ll|) ) compatible 
with the boundary conditions. The problem on the ring was considered in Ref. . 
Here one has asked which relations ( |1 . 1 1| ) are compatible with the trace operation 
and one has given the smallest representations (this was enough for applications). 

In Appendix A we present an instructive different (although equivalent) approach 
to obtain part of the diffusion algebras for A^ = 3. We also show a natural way to 
define a quotient of one of the algebras. 

In Sec. 4 we consider the case of A^ generators. First we give seven series of 
algebras and then present a theorem which allows to find all the diffusion algebras 
of PBW type. The proof of this theorem would imply a long discussion and would not 
fit in this paper which is also aimed at physicists who are not necessarily interested 
in combinatorics. We have therefore decided to publish it separately. Some algebras 
are not suitable for applications to stochastic processes (they are not compatible 
with positive rates). 

In order to help the reader who is interested in applications and not in mathe- 
matics, in Appendix B we list the diffusion algebras with positive rates for A^ = 4. 

In Sec. 5 several physically meaningful generalizations of the diffusion algebras 
are discussed. A possible connection between diffusion algebras and quantum Lie 
algebras is also pointed out. 



2 Diamond conditions. 

We shall search for PBW-type associative algebras which are generated by the unit 
e and the elements Pq,, a = 0, 1, . . . , A^ — 1 satisfying A^(A^ — l)/2 quadratic-linear 
relations given by (|1 . 1 1|) . 



The PBW property (see, e.g.. Sec. 3 of the Ref.[0]) implies that given a set of 
generators {"Po}, one can express any element of the algebra as a linear combination 
of ordered monomials in Furthermore, all the ordered monomials are assumed 
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to be linearly independent .0 

For concreteness let us fix an alphabetic order 

Vp>V^ if P>a. (2.1) 

Then, a linear basis for the PBW-type algebra is given by the unit e and the set of 
monomials 

Vl\VZ-...-Vll , k=l,2,... , (2.2) 

where ai > 02 > • • • > «fc and ni, n2, . . . , are arbitrary positive integers. 

Imposing the PBW condition for the diffusion algebra (|1.11|) we first demand 

ga^^O, V a < /? (2.3) 

in order to be able to express any polynomial in as a linear combination of the 
basic monomials ( |2.2|) . 

Next, using ( |1 . 1 1| ) one can reorder any cubic monomial 'D^'DpD-y —>■ 'D^'DpDa-, 



where a < /5 < 7 in two different ways: 

Demanding the coincidence of the resulting expressions for 'D^'DpD^ in terms of 
ordered monomials one obtains the relation 

+ 3^7 gpa{^a(3-f + ^Pa)T^pT^cx + XaXp^g^a — g^y/B — ^af3^f) T^-y 

+ Xo,Xy{g^P - gpo.)Vp + XaXpig^a - gya + ^aPy)1^a = 0, (2.4) 

which results in the following six conditions for the gaps and the XaS 

Xa g-yfsi^ay " ^ap) = , (2.5) 

Xpgyai^lSy + Kp) = , (2.6) 

X^ gpai^a-r - ^I3y) = , (2.7) 

XpX^ (A/3^ + gaf3 - ga-f) = , (2.8) 

XaXy {gpa - fi'7/3) = , (2.9) 

XaXp{Aap + gpy - gay) = , V « < /3 < 7 . (2-10) 



^ The term PBW is due to Poincare, Birkhoff and Witt |jT^ who describe a Unear basis in a 
universal enveloping algebra. 
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In Eqs. ( ^.41 ) and (|2.5|)-( PTTI1| ) we have introduced the notation 

:= da/S — Qpa , ^af3^ '■= ^afi + + A^q • (2-11) 

One can show that the conditions (|2.5| )- (|2?T0| ) are necessary in order to avoid 
hnear dependences between ordered quadratic (or even first order) monomials in 
Pq, Dp and T)^. One can improve this resuh applying the diamond Lemma |jT^ to 
our concrete case. Namely, the algebra ( |1 . 1 1| ) possesses the PBW property iff the 
conditions (|2.5| )- (|2.1CI| ) are fulfilled. We will refer relations (|2.5|) - (|2.1CI|) as "diamond 
conditions" in what follows. 

In the next sections we are going to find the solutions of the diamond conditions 
first for the case = 3 and then in the general case. 



3 Classification of diffusion algebras with 3 gen- 
erators. 

The classification of diffusion algebras generated by an arbitrary number A^ of ele- 
ments proceeds as follows. 

First, one notes that any subset oik < N elements generates a subalgebra of 
( |1.11| ) which is again a diffusion algebra. So, it looks natural to begin with minimal 



size subalgebras produced by 3 generators (for the case of 2 generators one does 
not have any nontrivial diamond conditions to solve). Then, a closer inspection of 
the diamond relations ( ^.5| )-( PTTI1D shows that they would be fulfilled for the whole 



algebra ( 1.11 ) provided that they are satisfied for all the minimal size subalgebras. 

So, we shall start by classifying diffusion algebras generated by three elements 
Pq, X>/3, X>^, a < /5 < 7. It is natural to fix the values of indices as a = 0, /? = 1, 
7 = 2 (this choice is adopted in Appendix A). Here however we will not assign 
concrete values to the indices a, (3 and 7 keeping in mind that for general A^ > 3 
a < /3 < 7 may denote any triple of indices from the set 0, 1, . . . , A^ — 1. 

Depending on how many of the parameters Xq,,X/3,x^ take nonzero values the 
classification falls into four cases. Namely, when all three parameters XaiXp and 

are nonzero we obtain the algebras of type A. If one of the x's is zero and the 
remaining two are not equal to zero we have the algebras of type B. The algebras 
with only one nonzero parameter x and those with all x's zero are called algebras 
of type C and respectively D. 

Now we consider in details the algebras of type A - D. 

Case A. All Xa^xp and x^ are nonzero. 
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Eqs.(P7^)-( |2.1(J| ) give us two constraints 

gpa = g-yf3 = Qafi + Qp^ — Qa-y ■ (3.1) 

Then, there are two possibihties. 

1). If g^oi = g-y/s 7^ 0, the Eqs.( p.5| ) and ( p.7| ) give A^^ = A^^ = A^^ whereof one 
obtains 

g-ya = g^-y = gaff ■ (3.2) 

In view of ( p.3|) one has g^a 7^ and, therefore, from eq.( p.6|) one obtains 

A/3^ = — Aq/3 . (3.3) 
Finally, from (p.!]), ( |3.2|) and ( p.3[ ) one concludes 

^ij = ^7^0, Vi, j G {a,/?,7} . (3.4) 
The corresponding y4/-type algebra is 



Ai 



g Pi^ ^j] = T^i - ' V i 7^ j G {a, /?, 7} , 5- 7^ . 



These are relations of Lie algebraic type. By rescaling the generators 
one can remove all the parameters from ( |3.5| ) and obtain: 

[Ei, Ej 



(3.5) 



^Ei 



Ei — Ej . 



2). If g^a = g-yi3 = 0, using Eq.( p.l|) one transforms the only remaining nontrivial 
equation ( p.6| ) as 

g-yai-^ISf + ^ap) = g^aig/Sy + gajs) = g-ya ga-y = ^ g.ya = . 

The corresponding A/^-type algebra is: 



A 



II 



gijViVj = Xj Vi-XiVj , \/ i < j e {a, P, 7} 
where gij := gi - gj , ^ gj y i ^ j . 



(3.6) 



The parameters gi {i G 7}) introduced in Eq. (|3.6| ) above are defined up to a 

common shift gi —>■ gi + c. 

This algebra is invariant under the transformation 

1 



Xi 



Xi 



gi-y 



. gi 



y-gi 



{y - giY ' 
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where y is an arbitrary parameter. As explained in Appendix A, the algebra 
has a natural quotient [P^, V^] = 0. We would like to mention that this algebra is 



already known ||T3|. 

Note that the algebras Aj and Ajj can directly be extended to the case i,j = 
0, 1, 2 . . . , — 1 for A^ > 3 (they correspond to the algebras Ai(N) and Ajj{N) 
discussed in Sec. 4). 

Case B. Among the coefficients Xa, xp and x^ one is equal to zero. 

Let Xa,x^ 7^ , X/3 = 0. In this case, the Eqs.(PTB|), (p^.8|), ( |2.1U| ) become trivial 
and Eq. (|2.9|) gives 

9f3a = ■ (3.7) 

There are two ways to satisfy the remaining Eqs. (|2.5|) and (|2.7|) . 

1). Eqs. (|2.5|) and ( p.7|) are satisfied if one chooses A^^ = A^^ = A^^ =: A which, 
in view of ( |3.7| ), leads to 



gaf3 — gp-y , g-ya — ga-y + g/Ba — gaf3 ■ 

The corresponding algebra is 



(3.8) 





gi3 VaVfS 


- {gp 


- A) VpVa = -XaVp , 


B(i) 


g VaVy - 


-{g- 






gp VpV^ 


- (gp 


- A) V^Vp = x^Vp , ^g, gp^Q . 



(3.9) 



Here for sake of future convenience we have parameterized the bulk rates via gf^, g 
and A. The algebra ( |3.9| ) is also known ||^, |T^. 

2). If g^ji = gpa = 0, then Eqs. (|2.5| ), (|2.7| ) are trivially satisfied and the algebra 
reads 





gal3 T^cJ^P 


= -x^Vp 




B(2) 




g^a T-^'yT-^a 






g(3^ VpV^ 


= x^Vp , 


V gal3, gay, gi3-i 7^ . 



(3.10) 



This algebra can be found already in Refs.[^, |T^. Notice that if one takes g^ = A 
in Eq. (p.9|) one obtains a special case of the algebra E^^) (Eq. (|3l0|) ). 

We will not repeat the same considerations for the 7^ 0, = and 

Xf3,x^ y^O,Xa = 0. The resulting algebras are: 
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' g V^Vp - ((? - A) VpV^ 

VaV^ - {g^ - A) V^Va 
. (5-7 - A) VpV^ - g^ V^Vp 

(ga - A) VaVfS - ga VpV^ 
g^ VJ)^ - {g^ - A) V^V^ 
^ g V(,V^ -{g-A) V^Vp 



-xpV. 



(3.11) 



7 ! 



Wgj^O, g,^{0,A} 



(3.12) 



,Vp-Xf3V^ , yg^O, ga^{0,A} ; 



B(3) 



g V^Vp -{g-A) V^V^, = 
{g^ - A) VpV^ = -xpV^ , 



XpVa - XcVp 



V(7^0,^7, ^{0,A}; 



(3.13) 



B(4) 



{g^ - A) V^Vf, 

g VpV^ -{g-A) V^Vp = x^Vp - xpV^ , V ^ , ^ {0, A} . 



X^T^a. 1 



(3.14) 



The algebras ( 3.11 ) and ( 3.12 ) are just different presentations of the algebra i?*^^) 
(one has to make the substitution /5 <-> 7 in Eq. (|3.11| ) and a <-> /5 in Eq. (|3.12| )). 

The relation between the algebras B^^^ and i?^'*-' is less trivial (therefore we keep 
them as different cases in classification). One can obtain the relations for the 5*^^-* 
algebra by inverting the order of all products in the i?*-^-* algebra (Eqs.( ^34|) ), i.e. 
by reading the relations ( ^.14| ) from the right to the left and changing the signs of all 
x's. That means, the algebras -B*^^-' and B'^^^ describe mirror (left-right) symmetric 
physical processes. 

The algebras 5 and 5^^) are completely different from B'^'^\ Even the number 
of independent rates in cases 5(3) and is not the same as in B^ K 

Case C. Two of the coefficients x^, Xf^ and x^ are equal to zero. 



In this case only one of the Eqs.( p.5| )-( p.lO| ) remains nontrivial and the analysis 
becomes straightforward. Below we present the relations for the type C algebras 



with Xa 7^ 0, Xfi 



0. The expressions for the cases xp 7^ 0, Xc 



= and 



7^ 0, Xq, = X/3 = can be obtained by the substitutions a 
a^7^/5^ain (|3J[5|) and ( gJGl) . 



[3 and, respectively. 
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9/3 T^J^p 


- {gp - A) VpV^ = 




C(i) 




- {g^ - A) V^V^ = 








- g-/i3 'D^'Dp = , 


gf3,g-y,gf3'y 7^ ; 



(3.15) 



C(2) 



gal3 T^aT^P — gpa T^pT^a 
gay T^aT-^y g^a T^yT^a 

gp^ VfsV^ = , 



XaT-^y 1 

V ga/3, gay, gi3y ^ 



(3.16) 



We observe that for g^f^ = 0, the algebra ( |3.15|) is a special case of the algebra 
( |3.16|) . For the sake of convenience (see Sec. 4), we will not stress this observation 
any further. 

Notice that in case A 7^ shifting the generator 
Eq. (|3.15|) one obtains 



Va — Xa/A, in 



' gp VJ)p - {gp - A) VpV^ = 
g^ VJ)^ - {g^ - A) V^V^ = 

. gf3y T^pT^y — gyP T^yT^P = , 



(3.17) 



which brings the C*-^-* algebra to the subcase of a family of quantum hyperplanes 
(see case D below)^ If one thinks of applications to stochastic processes, the shift 
we made is not an innocent one since for g/^^ = g^jj, the algebra (p.l5|) has represen- 
tations with traces (for example, the one- dimensional representation given by the 
shift) whereas the algebra ( |3.17| ) has none. 

A useful algebra (not of PBW type) is obtained if one takes not only g^yp = 
but also gi3y = in the algebra C*-^-* given by Eg. ( p. 161) . In this way one obtains the 
algebra used in Ref . [Q] . 

Case D. All the coefficients Xa, xp and x^ are equal to zero. 



In this case we obtain the algebra of Manin's quantum hyperplane [16 



sponding to multiparametric Drinfeld-Jimbo R-matrix 



D 



gab VaVb - gbaVaPb = , V a, 6 G {«, /3, 7} : a<h , gab ^ . 



corre- 



(3.18) 



Conversely, using the linear shifts of generators = Vi + Ui, \f i E {a,/?, 7} in the I?-type 



algebra (3.18) and demanding the resulting relations to agree with the diffusion algebra Ansatz 



( f .11 ) one recovers the C'^^^-type algebras only. 
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A different point of view in understanding some of the algebras presented liere 
is discussed in Appendix A. 



4 Diffusion algebras with > 3 generators. 



While classifying the PBW-type algebras ( 1.11 ) with more then 3 generators one 
meets the combinatorial problem of consistently combining several minimal subalge- 
bras generated by triples {VajV/^,!)^} (each of these subalgebras belonging to one 
of the types A-D listed in Sec. 3) to a larger algebra. In this section we shall first 
construct several basic series of diffusion algebras being extensions of the Aj, An, 
Q^T^f^ D -type algebras from the previous section. The C^^-'-type triples 
will appear later on in our considerations. There is a deep reason behind our choice 
of starting first with the algebra C^^-* and taking into account the algebra C^^-* later. 
In this way one can easier state the theorem presented at the end of this section 
and which is the central part of our work. A proof of this theorem will be given 
elsewhere. 

In Sec. 3 we have shown that the classification of diffusion algebras with = 3 
generators depends essentially on the number of nonvanishing parameters Xa in the 
Ansatz (|1.11| ). The same is true for general N . Therefore we shall split the set {a} 
labeling different species of particles (= different generators P^) into two subsets 
{q;}={z} U {a}. From now on we assign letters i, j, fc, etc. to the indices of the 
first subset and assume that Xj, Xj, x^, . . . 7^ 0. The indices of the second subset are 
denoted by letters a, 6, c, etc. and it is implied that = 0. Let 

A^i and Nq denote the number of elements of the first and second subsets. Clearly, 
A'o + A^i = — the total number of indices of both kinds. 

We should stress however that Ni — the number of nonzero x's — is the most 
noticeable but not the only relevant information for the classification. A supple- 
mentary information is given by a number of nonzero bulk rates Qap in the defining 
relations of the algebras (cf. cases Aj and Ajj, B^^^ and B^'^\ or C^^-* and C^^-* from 
Sec. 3) and the mutual arrangement of the indices {i} and {a} in the alphabetic order 
(see Eq.( p.l|) and the definition of the algebras B^'^\ B^^^ and B^^^ from Sec. 3). 

Algebras of type A. We shall start by considering the algebras with the number 
A^i of nonzero x's not less than 3 — we call them algebras of type A. Obviously, any 
such algebra contains a minimal subalgebra of type Aj, or Ajj. So, these algebras 
are naturally obtained by a sequence of consistent (in a sense of diamond conditions) 
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extensions starting with the = A^^i = 3 algebras ( |3.5| ) or ( |3.(j| ) and adding one new 
generator "Dq, at each step of the iteration. 

It is suitable to begin the extension procedure with the generators whose indices 
lie in the subset {i}. At the first step one adds a fourth generator, say T>i, to the 
triple {I)i,Vj,Vk} (recall, once again, that Xi,Xj,Xk,xi 7^ 0). The resulting algebra 
contains four minimal subalgebras of the types either Aj, or Ajj. An easy check 
shows that it is possible to combine only triples of the same type. Continuing the 
extension procedure one finally obtains algebras with A^^i generators for which all the 
minimal subalgebras are of the same type, either Aj or Ajj. The defining relations 
for these two types of diffusion algebras — Aj{Ni) and Ajj{Ni) — are given by the 
Eqs.(p.5|) and, respectively (p.6|), with i < j spanning the whole set {i}. 

We continue the extension procedure adding to the algebras Aiji{Ni), Nq new 
generators with their indices lying in the subset {a}. First, we add one generator, 
say Va, and take care that all the newly appeared triples {Va, T>i, Vj} belong to one 
of the algebras of type B (see Eqs.(|]9D, ( glOD , ( gl3D and (^A^). Next, adding a 
second generator, say V^, we again require that all the triples {I)h,Vi,Vj} are of 
type B and, moreover, demand that the triples {Va,'Di„Vi} are the generators of 
a C/-type algebra ( p.l5|) . Adding new generators we have also to impose D-tjpe 
algebraic relations for the triples {Va,T>h,Vc}. As a result, we obtain two different 
extensions for the algebra Aj{Ni) 





9 Vj] -- 




e {i} : j- 


97^0, 




ga [Va, Vi] 




Wi e {i}, Va G {a}; 






gab VaVh - 


QboPhVa = , 


^ a,b & {a} : a < b; 


Qab^O ■ 



(4.1) 



9 [V,, V,] 


= Xj Vi - 


x,Vj , y i,j e {i} : i ^ j; 


97^0, 


g+ Vi Vb 


= -Xi Vb , 


V(i G b G {a}) : t 


<b; g+^ 0, 


9- VaVi 


Xj Vfi , 


V(i G {i}, a G {a}) : i > a; g_ , 


where g. 


^ = -9+ , 


if there exist Di, Da, Db : 


a < i <b , 


9ab VaVb 


- gbaVbVa 


= , W a,b E {a} : a < b; 


9ab1^^ ■ 



(4.2) 



Here and in what follows we indicate two integers A^o and A^^i — the numbers of 
indices in the sets {i} and {a}, respectively — in braces to specify the type of 
algebra. As the reader can notice looking closely at Eq. ( |4.2| ), specifying Nq and 
Ni one obtains several algebras all denoted by A\ {Ni, Nq). We didn't introduce a 
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different notation for eacli algebra in order to simplify the notations. We adopted 
the same attitude also for other algebras described below (see Eqs.(^73|) and (^.5|)). 

The algebra Afj{Ni) possesses a unique extension 



Aii(iVi,iVo 



Qij j '^i J 


y i,j e {i} : i< j 




gi+ T^i T^b = -Xi Vf, , 


V(2 G {i}, b G {a}) 


: i <b, , 


Qi^ Da Vi = XiVa , 


W{i G {i}, a G {a}) 


: i > a . 


Here gij ■= g^- , gi^ 


\--=g+ + gi, gi- ■= 


g--gi , 


where for all i: gi^ —g^ 


, giT^ g- ; giT^ gj 




and g- = —g^ , if there exists Di,Da,Db : 


a < i < b . 


gab T^J^b - gbaPbT^a = , 


\f a,b e {a} : a < b; 


gab^o ■ 



(4.3) 



The parameters g^ and g^ in algebras Af'\Ni,NQ) and Aji{Ni, Nq) remain 
independent provided that mutual order of indices of the subsets {i} and {a} (and 
hence the order of the generators Di and Da, see (p.l|)) is like follows 



ii<i2< ... <ik < ai < a2 < . . . < a^o < ik+i < ik+2 < . . . < ini ■ 

Only in this case all the 5-type minimal subalgebras in the algebras ( [4.2|) and (^4.3|) 
belong to the type B^'^\ In the presence of B^^\ or B^^^ type triples the parameters 
and g^ are constrained by condition g^ + g^ = (keep in mind that for stochastic 
processes all rates have to be non-negative). 

Algebras of type B. We now consider the case where the set {i} contains exactly 
two indices (A^^i = 2), say, i and j, i < j. To obtain such algebras — we call them 
the algebras of type B — one should consider consistent extensions of the i?-type 
triples ( |3.9| ), ( |3.10| ), ( p.l3|) and ( p.l4| ) by {Nq — 1) generators with their labels in the 
set {a}. 

Starting with the .B^^^-type triple {Vi,Vaj^,T)j}, the only possibility is to add 
new generators Va^, T^a-n etc. such that all the minimal subalgebras {Vi,Va2,T^j}, 
{Vi,Va._i,Vj}, etc. satisfy again i?*^^^-type relations (|3J| ). In this situation, we can 
arrange the alphabetic order of the generators as follows: Vj > Va^ > Va^ > . . . > 



T^aj^^ > 'Di (cf. with the comment below Eq. ( |3.14D ). So it is natural to put in this 
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case i = 0, j = N — 1. The extended algebra B*^-'-)(2, A^^o) reads 



g Do Vn- 


I - (g- A) Vn-1 Vq = xn-1 


Vq — Xq Vn-1 5 


97^0, 


ga Vq Va 


- (ga - A) Va Vq = -Xq Va , 






ga Va Vn 


-1 - (ga - A) Vn-1 Va = XN- 


1 Va , Vl<a< AT-l, 




9ab Va Vb 


- gba VbVa = , 


W a,b E {a} : a < b; 


9ab^0 ■ 



(4.4) 



The situation becomes different for the other 5-type triples ( |3.10D , ( p. 131) and 



.141). Extending these algebras one can get algebras containing all the B^'^\ B^^^ 



and i?'^^)-type minimal subalgebras. Therefore we introduce a unified notation 
B^'^\2, «<, n, n>) for extensions of the triples B^'^^ — B^^\ Here n< + n + = Aq 
and a mutual order of the indices i, j and the indices from the set {a} is as follows 

Oi < . . . < < i < a(n<+l) < . . . < a(n<+n) < j < «(n<+n+l) < • • • < UNq ■ 

The algebra B^^\2, n<, n, riy) reads 



gViVj - {g - A) V, Vi = 




97^0, 


g+ Vi Va = -Xi Va , 


{g+ - A) Vj Va = -Xj 


Va, \f a>j , 


ig--A)VaVi = XiVa , 


g- Va Vj = Xj Va , 


W a < i , 


g+ Vi Va = ~Xi Va , 


g- Va Vj = Xj Va , 


V i < a < j , 


/ if n< 
where g+ t S , 

1 A if ?2> 


< A^o / / if 
>0 ' ^-^1 A if 


n> < No 
n< > 


and (7+ + (7_ = A if among 


the numbers n, riy there are two nonzeros, 


gab VaVb - gbaVbVa = , 


\/ a,b E {a} : 


a<b; gab 7^0 . 



(4.5) 



The algebra B^'^\2, n^,n,ny) contains i^'^^^-type minimal subalgebras ({Va,Vi, 
Vj} for a < i), n B^'^^-tj'pe minimal subalgebras {{Vi,Va,Vj} for i < a < j) and 
n> 5(3)-type triples {{Vi,Vj,Va} for a > j). 

Algebras of type C. Next, we consider algebras with only one index i in the subset 
{i}, let us call them the algebras of the type C. These algebras arise from extension 
of the Ci triple {V^, Va, Vb} by (A'o — 2) generators V^, . . . labeled by indices from 
the subset {a}. Checking the consistency of such an extension is straightforward 
and, therefore, we shall just present directly the resulting algebra 
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- {ga - A) Va Vi 




C(l,iVo) 


where 


ga^O if i <a , 


and if a < i , 




gab 'DaVb 


- gbJ^bDa = , 


\/ a,b e {a} : a < b; gab 7^ ■ 



(4.6) 



As in case of the C(i)-type triple (|3J[5|) , for A ^ one can reduce the algebra (|4.6|) 
to a subcase of the family of quantum hyperplanes (see case D below) shifting the 
generator X^j: Vi T>i — Xij K. Nevertheless, as will be shown in the theorem given 
below, it is useful to keep the definition given by Eq. (|4.6| ) for the algebra C(l, A^^q) 
since in this way one can use it as a building block for the construction of new 
algebras. In the new algebras the shift will not be possible anymore. 

Note that, unlike all the previous cases, one can consider the algebra C(l,l) 
produced by a pair of generators. We shall use this possibility in the theorem stated 
below. For instance, the C^^-* algebra (|3.16| ) can be constructed as a combination of 
two C(l, 1) algebras by a blending procedure described in the theorem given below. 
Further examples of an application of C(l, 1) algebras are given in Appendix B (see 
the cases 17 and 18). 



Case D. The D algebras — the algebras with no indices in the subset {i} 
represented by a family of A^ = A^o dimensional quantum hyperplanes 



D(0,Aro 



ga 



^VaVi, 



gba'DbVa = , \f a,h e {a] : a < h] gab ^ . 



are 



(4.7) 



Now we are ready to complete a classification scheme. To do so one needs to take 
into consideration the possibility of using C^^-'-type triples in the algebra extension 
process. As a result one derives a procedure to obtain all the the diffusion algebras 
which is described in the following theorem. 

Theorem. A diffusion algebra has Ni generators Vi, where one assumes Xi ^ 0, 
(i = 1,2,..., Ni) and Nq generators Va, with = 0, fa = 1, 2, . . . , Nq). If we 
don't distinguish between the two kinds of generators, we denote them by Va, (a = 
1,2,...,N = No + N,). 

IfNi = 0, the algebras are D{0,No) (see Eq.(^)). 

If Ni 7^ 0, all diffusion algebras can be obtained by a blending procedure using 
the algebras l \4-l\ )- i^4-A )- The blending procedure can be described as follows. 

Consider two of the algebras ^)-^), denoted byX{Ni, N^''^) andY{Ni, N^^^), 
both having the same number of generators Vi which satisfy the same relations among 
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themselves in the two algebras, and generators T>a^ (respectively T^ay)- Through 
blending, one can obtain a new diffusion algebra Z{Ni,N^''^ + N^^^ = No) with 
generators V, and Va, ("a = 1, 2, ... , Nq). Since the X{Ni, N^""^) and Y{Ni, N^^^) 
algebras are both of PBW type, the N^^^ indices (respectively the Nq^^ indices) are in 
a given alphabetic order. We blend now the A^q^'* and Nq^^ indices together in an ar- 
bitrary alphabetic order but respecting the order for the Njf^ indices ( respectively the 
Nq^^ indices) which are fixed in the algebra X{Ni, N^^^) (respectively Y{Ni, Nlf'^)). 

(x) 

The alphabetic order of the Nq indices in respect to the Ni indices (respectively the 
A^o^^ indices in respect to the Ni indices) given again by the two algebras X{Ni, Njf^) 
and Y{Ni, Nlf^), has also to be respected. For each alphabetic order of the Nq indices 
one obtains a new algebra in the following way. The relations among the genera- 
tors Vi and those among the generators T>i and (respectively among T>i and 
T^ay) coincide with the relations in the algebras X{Ni, Nq^^) and Y{Ni, N^^^) . The 
remaining relations among the generators T>a^ and T>ay are: 

l^a.'Day = , \f a^,ay: < ay , 

'Day'Da^ = , 'i a^.ay-. > ay . 



A "blended" algebra can now be blended with one of the algebras (^^-(^7^ and one 
can obtain a new algebra. 

It is important to stress that through the blending procedure one can obtain the 
same algebra using different blendings. Therefore what we have is a construction 
rather than a classification of the diffusion algebras. 

We would like to point out that it is easy to see that one can blend together only 
algebras of the same type: A^^ or Af^ with algebras A'-p or Af^ , An with An, B^^^ 
or 5^^) with B^^^ or B'^'^\ and C with C. One can show that blending together two 
A^P (respectively 5*^^^) algebras leads to A^P (respectively S*^^)) algebras. Therefore 
one can use an A^ (respectively 5*^^^) algebra only once during blending procedure. 
At the same time the algebras Af^ (respectively B^'^^) can be blended any number 
of times. 

Let us show an application of this theorem. Consider the algebra 

QQa VoVa - gaoVaVo = -XoVa , a = 1, 2, . . . , A^ - 1 , (4.9) 
QabV.Vh = , V a, 6 : 1 < a < 6 < AT - 1 . (4.10) 

This algebra is obtained taking A^ — 1 copies of the algebra C(l,l) ({«} = {0}, 
{a} = {1, 2, .., A^— 1}) and blending them together. Taking A^ = 3 and 4 one recovers 
the algebras given by Eq. (|3.16|) respectively Eqs.( p.25|) , ( p.26| ). If in Eq.( |4.10| ) one 
takes the rates gab = one recovers the algebra discussed in Ref.[@]. 
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5 Discussion. 



We have defined diffusion algebras. Those are PBW algebras with generators 
satisfying the relations ( |1.11|) . These algebras are useful to find stationary states of 
the stochastic processes given by the rates Qa/s- For the N species problem one finds 
several series of algebras which might be useful in applications using the Eq.( |1.3| ) 
for a ring and the Eqs.( |1.7|) and ( |1.8|) for open systems. Much work is still left. For 
example one has to find which boundary conditions, if any, are compatible with each 
of the algebras. We would like to stress that all the cases which were used up to 
now for applications are special cases of our construction. 

An open and relevant question is: do non PBW type "physically meaningful" 
algebras satisfying Eq.( |l.ll| ) exist .0 Such algebras could eventually represent "ex- 
ceptional algebras" similar to those which appear in the theory of simple Lie algebras 
or super algebras. 

An interesting and unsolved problem is the connection between the diffusion 
algebras and the so called quantum Lie algebras (see e.g. ||l8l) which look similar 
and which are not fully investigated and classified. In order to show that such a 
connection might be possible we will give here two examples of quantum Lie algebras 
related to the so-called Cremmer-Gervais R matrix [OT 



gp VaVp - {gp - A) VfiV^ = -x^Vp , 

gVJ)^-{g- k)V^V^ + w iVpY = XyVa - XaV^ , {5.1] 

gp VpV^ -{gp- K) V^Vp = x^Vp , 



and 



g^^ VaV.^ - g^a + w {VpY = x^V^ - xJD^ , (5.2) 

gp^ VpD^ = x^Vp , 

They are clearly extensions of the 5*^^^ and 5*^^^ algebras to which the term {T^pY 
has been added {w is an arbitrary parameter). The algebras (|5.1| ) and ( |5.2| ) are 
of PBW type. Quantum Lie algebras can also be relevant in a different context 
un- related to stochastic processes. As pointed out in Ref . pO| , quadratic algebras 



are useful also to describe the ground-states of one-dimensional quantum chains in 
equilibrium statistical physics if it happens that the ground-states have energy zero. 
This is a whole which is worth exploring. 

Before closing this paper, we would like to mention a natural extension of 
our results. The starting point of our investigation were the processes given in 



These algebras can be finite or infinite dimensional. 
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Eq.( |l.l| ) which are related to quantum hyperplanes (take all the x's equal to zero 
in Eq.( |l.ll| )). One can consider more general stochastic processes in which the bulk 



rates are related to quantum superplanes [|^. In these cases one obtains equations 



which generalize Eq.( l.ll|) (see Ref.|P]). This would lead us to something to which 



one could coin the name of "reaction-diffusion algebras". For the time being this is 
not more than a nice thought. 
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Appendix A. Comments on diffusion algebras with 
three generators. 

We consider the relations (|1.11| ) in the case when we have only three generators: 



920^2^0 - 5'02'Z^0^^2 = Xo "^2 - ^2 Vq , (A. 1) 

gi2'DiV2 - g2iT>2'Di = X2 "Ci - xi V2 , 
and look for the case when the rates satisfy the condition: 

A012 = goi - gio + gi2 - g2i + g2o - go2 = (A.2) 

It is trivial to verify [12] that the 1x1 matrices (c- numbers) 

A = y, (^ = 0,1,2), (A.3) 

Ji 

where 

/i = /o + 5-01 - gw , /2 = /o + 5-02 - g2o , (A.4) 



verify the relations ( |A.1|) (/o is an arbitrary parameter). One can use ( |A.3D in order 
to compute, using Eq.(|1.3|), the probability distribution on a ring. The probability 
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distribution one obtains for the stationary state is trivial (one has no correlations). 
The physics of the stationary state can however be interesting if one takes an open 
system. We then have to use Eq. (|1.8| ) and the one-dimensional representation of 
(|A.1|) is not of much help. If the rates satisfy only the condition ( A. 2 ) it is not clear 
if one doesn't have only the one-dimensional representation. At this point one can 
understand why we are interested in algebras of PBW type. For algebras of PBW 
type, we can be sure that one gets other representations (at least the regular one). 
The price to pay is that we will get more constraints on the rates than those given 
by Eq. (|A.2| ). In order to get algebras of PBW type compatible with the relation 
( [A. 21 ), it is useful to write the diamond condition ( ^.4| ) in a different way: 



(A.5) 



A012 {xq g2i D2 Di + xi g2o D2 Dq + X2 921 D2 Di + Xi X2 Dq - xq Xi D2) 



+ xogi2g2i [-Di, D2] + xigQ2g2Q [-D2, -Do] + a;2fi'oi5'io [-Do, D^i] = 
We now take into account Eq.( [A.2| ). There are several solutions of the Eq.( [A.5D 



a) gij = g . 

This gives (using Eq.( |A.l|) ) the algebra Aj (see Eq(3.5)). 



(A.6) 



b) 




920 = 9io = 
901 + 912 ■ 



(A.7) 



This gives the algebra Ajj (see Eq.( |3.6|) ) 



920 = 910 = , 



(A.8) 



This gives actually again the algebra Ajj, defined in Eq. (p.6|) with the substitution 

6-12 - 6-21 — ' gi2 ■ (A.9) 

This different derivation of the algebra Aji has a bonus: we have learned that we 
can take the quotient given by the second relation in Eq.( |A.8p of this algebra. 



d) 



Xo = , 
901 = 902 = , 
9i2 = 9 , 921 = 9- , 
. fl'io = 5'7 ' 920 = 9-1 - ^ 



(A.IO) 
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which is the algebra B^^^ (see Eq.( ^.13D ). 



xo -- 
910 

9l2 

901 



0, 

= 5-20 = 
-- 9 , 921 

-- 9a- : 



902 = 9a , 



(A.ll) 



which is the algebra B^^^ (defined in Eg. ( p. 141) ). 

Since the algebras Aj, An, B^^^ and 5*-^^ are all derived starting from Eq.( |A.l|) 
with the conditions ( |A.2| ), all these algebras have at least a one-dimensional repre- 
sentation. When the classification of the relations ( [A.l|) for which a trace operation 



exists was done |12[, there was no need to consider them separately since for the 
trace operation it is enough to have the expressions (|A.3[) and (|A.4|). 



Appendix B. Diffusion algebras with positive rates 
for iV = 4. 

Here we list all the diffusion algebras with = 4 generators which can be useful 
for stochastic processes with four species of particles (one can choose all the rates 
non-negative). 

We first describe those of the algebras from the seven series described in Sec. 4 
(see Eqs.( [4.1| )-( [4.7| )). For the reader's convenience we point out all the = 3 
subalgebras for each algebra in the list. The reader may be surprised by the fact 
that in the list which follows the same symbol will denote two algebras (see Eqs. 
( p.3| ) and ( |B.5| )), this is due to the fact that as can be seen already in Eqs.( |4.2| ) and 
([4.3|) the same notation is used for several algebras. 



0,1,2,3. 



:b.ii 



1. A/(4). 

g [Di, Vj] = Xj Vi - Xi Vj , i,j 

All the A^ = 3 subalgebras of this algebra (i.e., the subalgebras generated by 
different triples of the generators Vi, i = 0, 1, 2, 3) are of type Aj (|3.5|). 

2. A«(3,l). 



g [Di, Vj] = Xj Vi - Xi Vj , i,j = 1, 2, 3 
go [Vq, Vi] = XiVo . 



:b.2) 
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This algebra contains the Aj subalgebra {T>i^T>2-,T>^} and three B^^^ subalge- 
bras ( |3.9| ) in which one takes A = and which are generated by the triples 
{Vq,Vi,V2}, {Vo,Vi,Vs} and {Vo,V2,Vi}. 

l(2) 



3. A} '{3, 1) (two algebras). 



g [Vi, Vj] = Xj Vi 
go Vq Vi = XiVo , 

I, J = 1,2,3. 



and 




-xi Do , 

-- X2,3 T>o , 



(B.3) 



1,2,3. 



Both algebras presented in ( [B.3|) contain the Aj triple {T>i,T>2,T>^} and the 
triple {Vq,V2,V^} (see (|37[4l ) specialized to A = 0). The triples {Vq,Vi, 
V2] and {VQ,Vi,V-i} are of -B^^^ type for the first case in ( |B.3| ) and they 
belong to the type i?*^^^ for the last case (see (p.lOp with the special choice 

ga-y g^a fi')- 

There are two more algebras of type A\ [3,1). Their defining relations are 
obtained by inverting the order of the generators in all the products in formu- 
lae (p.3| ). This means that the last pair of algebras would describe physical 
processes which are just mirror reflections of the processes corresponding to 
the algebras 



4. An{A). 

gij T>i Vj = Xj T>i — Xi T>j , where i,j 



0,1,2,3 , and i < j . (B.4) 



Here gij := gi — gj and go > gi > g2 > gs so that gij > for all i < j. 
The algebra ( B.4 ) contains only Ajj subalgebras (|3.6|) . 



5. Ajj{3,l) (two algebras). 



gij Vi Vj 
goi Vo Vi = XiVo , 



and < 



gij Vi Vj 
g+ Vi Vq = -xi Vq , 
g^2 1^0 T^2 = X2V0 , 
5(_3 Vo V3 = X3V0 . 



(B.5) 



Here both i and j take the values 1, 2, or 3; furthermore one has z < j in the 
first relations in (p^ ); gij := gi - gj, goi := go - gi and g^i := g^ - gi, where 
9o> 9i> 92> 93 and g^ > 5(3. 
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The algebras in ( [B.5| ) contain the Aji triple {Di, 1)2,^^3} and the B^^'^ triple 
{Vo, V2, V3} (see ( CT ) with ^ = A). The triples {Vq, V^, V^} and {V^, V^, V3} 
are of B^^^ type for the first case in (p.5| ) and they belong to type -B*-^^ for the 
last case (see ( p.lO[) with g^a = 0). 

One can construct two more algebras which are mirror partners of those listed 
in era (c.f. the case 3). 



6. 5W(2,2) 



gVoVs - ((7 - A) I?3 Vo = X3 - xo V3 , 

Qa Vo Da - (ga " A) Da Vq = -Xq Va , 

ga Va P3 - {ga " A) V3 Va = X3 Va , tt 

gi2 Vi V2 - g21 V2V1 = . 



1,2, 



:b.6) 



This algebra contains two B^^^ triples {Vo,Vi,Vs} and {Vo,V2,Vs}, 

and two C(i) triples (|37[5|): {Dq, ^^i, ^^2} and {Pi,p2,^^3}- 



7. 5(2)(2, 0,2,0). 



(7 Do I?3 - (^7 - A) P3 ^^0 = 2:3 - a;o V^ 

g+ Vo Va = -Xo Va , 

g^ VaVi = xsVa , a = 1,2 , 

gi2 Vi V2 - g2i V2V1 = . 



(B.7) 



This algebra contains two E^^) triples (|3lOD : {Vo^Vi^V^} and {Pq, ^^2, ^^3}, 
and two C^^^ triples: {Vo,Vi,V2} and {1)1,1)2,^^3}, where in Eq.( p.l5| ) one 
takes g/3 = g-y = A, respectively gp = g-y = 0. 



8. 5(2)(2, 2,0,0). 



(? ©2 I?3 - (^7 - A) P3 ^^2 = 3^3 2^2 - a;2 V3 

{h - A) P„ V2 = X2Va, 

h Va V3 = X3 Va , a = 0, 1 , 

goi Vo Vi - gio Vi Vq = . 



(B.8) 



This algebra contains two B^^^ triples (^J§: {Pq, ^^2, ^^3} and {Vi,V2,V3}, 
and two C^^^ triples (|3l5| ) with g^ = g^ = 0: {©0,^^1,2^2} and {Di, ©2, ^^3}- 

The algebra 5(2) (2, 0,0,2) is a mirror partner of the algebra above. 



24 



(B.9) 



9. C{1,3). 

Qa Vo Va - {ga " A) Va Vq = -Xq T>a , 
gab Vb - gba l^b l^a = , O, 6 = 1, 2, 3 . 

This algebra contains three C^^^ triples {Vq,Vi,V2}, {Pqi^Ij^^s} and {V^, 
V2,V2}, and the triple {Pi,p2,2^3} of the type D ( glSl) . 

For A 7^ 0, one can do the shift Vq ^ Vq — xq/ A and bring this algebra to the 
subcase of the family of quantum hyperplanes ( [B.10|) . 

10. D(0,4) algebra, or quantum hyperplane. 

gabVaVb-gbaVbVa =0 , a, 6 = 0, 1, 2, 3 . (B.IO) 
Obviously, all the triples here are of D type. 

Next, we use the procedure of blending several diffusion algebras as described in 
the Theorem at the end of Sec. 4. There is no need to point out anymore the = 3 
subalgebras for each example separately, since the pair of main = 3 constituents 
which are blended to produce an = 4 algebra are mentioned explicitly in each 
case. The remaining two = 3 subalgebras always belong to the type C*-^-* ( |3.16| ). 

Note that blending algebras of the types Aj and An produces only examples 
with A^ > 5. 

One can get A^ = 4 diffusion algebras by blending any two of the following 
A^ = 3 type B algebras: B'-^\2, 1) {= B^^^ in the notations of Sec.3), 5(2)(2, 1, 0, 0) 
(= 5^), 5(2)(2, 0,1,0) {= 5(2)), and 5(2)(2, 0,0,1) {= B^^^). The results are 
listed below. 

11. Gluing B^^^ (with generators {T)q, Vi, V^}) and B^'^^ (with generators {Vq, D2, 
"Da}) one obtains 

gVoV; -{g- A) P3 = x-, Vo - Xo , 
g^ Vo V, - {g, - A) V, Vo = -Xq V, , 
g^ Pi P3 - {g^ - A) V3 Pi = S3 Pi , 

^ g+ Vo P2 = -Xq P2 . g- V2 P3 = X3 P2 , 

This set of relations should be supplemented by the condition 

gu V,V2 = . (B.12) 



The algebra (p.ll| ), (|B.12|) has a mirror partner with an opposite order of 



indices 1 and 2 (one uses the condition 5'2iP2^^i = instead of ( [B.12|) for it). 
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12. Gluing ( with generators {Vq, Vi, V^}) and 5^^^ (with generators {Do, ^^2, T^i}) 
one obtains 

' ^ Po 2^3 - - A) 2^3 ^X3Vo- xo Pa , 

< g+ Vo = -Xo , g_ V^, = , (B.13) 

^ h+ Do V2 = -Xo V2 , h_ V2 = X3V2 , 
and gi2 Vi V2 = . (B.14) 

13. Gluing (with generators {Pi, P2, Pa}) and S^^^ (with generators {Vo, V2, Pa}) 
one obtains 

^ P2 Pa - (^? - A) Pa P2 = xa P2 - X2 P3 , 
(/i - A) Po P2 = a;2 Po , /i Po Pa ^ X3 Vq , 
{gi - A) Pi P2 - ^1 P2 Pi = X2 Pi , 
gi Pi Pa - (^1 - A) Pa Pi = xa Pi , 



(B.15) 



and either goi Pq Pi = , or gio Pi Po =0 



(B.16) 



Gluing the algebras B^^^ and B^^^ produces a mirror partner of this algebra. 

14. Gluing B^^^ (with generators {Pi, P2, Pa}) and B^^^ (with generators {Po, P2, Pa}) 
one obtains 



' ^ P2 Pa -{g-A) Pa P2 = xa P2 - X2 Pa , 

{h - A) Po P2 = X2 Po , h Po Pa = xa Po , 
[ (/ - A) Pi P2 = X2 Pi , / Pi Pa = xa Pi , 



(B.17) 



and ^01 pQ Pi = . (B.18) 
Gluing the algebras B^^^ and B^^^ produces a mirror partner of this algebra. 

15. Gluing (^ith generators {Pi, P2, Pa}) and S^^) (with generators {Po, Pi, P2}) 
one obtains 



^ Pi P2 -{g- A) P2 Pi = X2 Pi - xi P2 , 

Pi P3 = -xi Pa , {h- A) P2 Pa = -X2 P3 , 
(/ - A) Po Pi = xi Po , / pQ P2 = a;2 pQ , 



and 510a Pq Pa =0 



(B.19) 
(B.20) 
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16. Gluing B^'^^ (with generators {Vi, T)2, ^^s}) and -B'-^-' (with generators {Vq, Vi, T>^ 
one obtains 

' gViV-i -ig-A) P3 Vi = Vi - xi , 

< g+ Pi I?2 = -xi V2 , (7- ©2 = ©2 , (B.21) 
^ (/i - A) Po 1^1 = Xi Po , 2^0 ^^3 = X3 , 

and ^02 ^^0 = . (B.22) 
Gluing the algebras i?*^^-* and B^^^ produces mirror partner of this algebra. 

Gluing a pair of -B*^^^ algebras gives a special case of the 5(^)(2,2) algebra see case 
6 above) with (721 = 0. 

There are two other possibilities to blend the C type algebras C(l, 2) and C(l, 1) 
(for their definition see Eq.( [4.6[ )) into a = 4 diffusion algebra. 

17. Gluing C(l,2) (with generators {Vq,Vi,V2}) and C(l,l) (with generators 
{Vo,!)^}) one obtains 

ga Vo Da - {ga - A) 'Da Vo = -xo Va , a = 1, 2 , 

< (712 I?2 - g2i V2V1 = , (B.23) 
. gm "Po ^^3 - 5-30 ^^3 ^^0 = -Xo "^3 , 

with either one of the following two sets of conditions: 

(713 VJ)-, = (723 V^V^ = , or (713 P1P3 = (732 V^V2 = . (B.24) 

In this algebra, besides the C^^^ triple {Vq,Vi,V2} there are two C^^^ triples 
{1)0,2^1,2^3} and {Vo,V2,V3} and the D triple {Pi, ©2, 2^3}- 

18. Gluing three copies of C(l, 1) algebra one obtains 

goa Vq Va - gao T^a 2^0 = "Xq 2)^ , a = 1, 2, 3 , (B.25) 

and (712 2)iP2 = ^7i3 2^12^3 = ^723 2^22^3 = . (B.26) 

This algebra contains one D type triple {Di, 1)2, 2^3} and all the other triples 
are of the type C^'^\ 
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